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Abstract 

The neutrinoless double /9-decay is reviewed. Model independent 
evidence in favor of neutrino masses and mixing is briefly summarized. 
The data of the recent experiments on the search for Oi^/3/3-decay are 
presented and some future experiments are discussed. The possible 
values of the effective Majorana mass, which can be predicted on the 
basis of the neutrino oscillation data under different assumptions on 
the pattern of the neutrino mass spectrum, are considered. A possible 
model independent test of the nuclear matrix element calculations is 
discussed. 



1 Introduction 

The status of the problem of the neutrino mixing drastically changed during 
the last 5-6 years: in the atmospheric Super-Kamiokande experiment PP, 
the solar SNO experiment j2l 0] and the reactor KamLAND experiment ji] 
strong model independent evidence of neutrino oscillations was obtained. 

There are many open problems of neutrino mixing. The value of the pa- 
rameter sin^ is the most urgent one. The value of this parameter is crucial 
for the search of such fundamental effects of the tree-neutrino mixing as CP 
violation in the lepton sector. The value of the parameter sin^ ^13 will be 
measured in the future reactor and long baseline accelerator experiments 

One of the most important problem of the neutrino mixing is the problem 
of the nature of the neutrinos with definite masses uf. are they truly neu- 
tral Majorana particles or Dirac particles, possessing conserved total lepton 
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number? The solution of this problem will have a profound impact on the 
understanding of the origin of small neutrino masses and neutrino mixing. 
The study of the neutrinoless double /?-decay (Oi//3/3-decay) 



is the most sensitive method of the investigation of the Majorana nature of 
neutrinos with definite mass. We will discuss here this process. In the next 
section model independent evidence of neutrino oscillations will be briefly 
reviewed. Then the data of the recent Oi//?/?-decay experiments and projects 
of the future experiments will be discussed. We will consider the possible 
values of the effective Majorana mass, which can be predicted from the results 
of the analysis of the neutrino oscillation data. In the last section we will 
discuss the problem of nuclear matrix elements of the Oz//3/3-decay. A possible 
test of the models of the nuclear matrix element calculations will be proposed. 

2 Model independent evidence of neutrino 
oscillations 

Compelling model independent evidence of neutrino masses and mixing was 
obtained in the atmospheric Super-Kamiokande P, solar SNO [21 El and 
reactor KamLAND |1] experiments. 

In the Super-Kamiokande (SK) experiment zenith angle 9z dependence 
of the electron and muon atmospheric neutrino events was studied in details. 
If there are no neutrino oscillations the number of the high-energy neutrino- 
induced /i(e) events N^^e{cos6z) must satisfy the symmetry relation 



For the high-energy muon neutrinos this relation is strongly violated: the 
significant deficit of the up-going muons is observed. For the ratio of the 
total numbers of the up-going and down-going muon neutrinos the following 
value 



was obtained. Here U is the total number of the up-going muons (—1 < 
cos^^ < —0.2) and D is the total number of the down-going muons (0.2 < 
cos ^2 < 1). 



{A,Z) + + e~ + e 



iV^,e(cOs6'^) = A^^,e(-COs6'2). 




(1) 
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The up-going muons are produced by neutrinos passing distances from 
about 500 km to about 13000 km and the down-going muons are produced by 
neutrinos travehng distances from about 20 km to about 500 km. The ratio 
(P) clearly demonstrates the dependence of the flux of the muon neutrinos 
on the distance between the neutrino production region in the atmosphere 
and neutrino detector. 

The data of the SK experiment and other atmospheric neutrino exper- 
iments MACRO Hi and SOUDAN 2 jH] are described by the two-neutrino 
t'/i oscillations. From the two-neutrino analysis of the SK data the 

following best-fit values of the oscillation parameters were found 

(Am2)sK = 2 ■ 10- (sin^ 2^)sk = 1.0 (xLn = 170.8/170 d.o.f .) (2) 

At the 90% CL the oscillation parameters are in the ranges 

1.3 • 10"^ < (Am2)sK < 3.0 ■ 10" (sin^ 2^)sk > 0.9 (3) 

In the SNO experiment [21 E] solar neutrinos from the decay *B -^^ Be -|- 
e"*" -|- z/g are detected via the observation of the CC reaction 

Ve + d^e^+p + p (4) 

and the NC reaction 

vi + d-^vi + n + p (/ = e, yU, r) (5) 
For the total flux of z/g detected via CC reaction Q the following value 

= (1.59^H?(stat.)^S:S^ (syst.)) ■ 10^ cm"^."^ (6) 

was obtained in the SNO experiment. 

For the total flux of z/g, z/^ and Vr detected via NC reaction (jSJ the sig- 
nificantly larger value 

$SNO _ ^5 21 ± 0.27 ± 0.38) ■ 10*^ cm'^s^^ (7) 

was found. Eq.® and Eq.((7j) give us a model independent evidence of the 
transitions of the original solar z/g into v^^t- 

The results of all solar neutrino experiments [3 El 13 UHl HH are 
described by the two-neutrino z/g z/g survival probability in matter. From 
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the analysis of the data in the most preferable LMA region the following 
best-fit values of the oscillation parameters were found [2] 

(Am'),oi = 5 ■ 10- (tan^ 9U = 0.34; xLn = 57/72 d.o.f .. (8) 

In the KamLAND experiment electron antineutrinos from many reac- 
tors in Japan and Korea are detected via observation of the process 

t'e + P — ^ e"^ + n 

The average distance between reactors and the detector in this experiment 
is 180 km. For the ratio of the total numbers of the observed and expected 
events the following value 

= 0.611 ± 0.085 ± 0.041 

^exp 

was obtained. This result is a clear evidence of the disappearance of i/g's on 
the way from the reactors to the detector. 

The KamLAND data are described by the two-neutrino z/g Uf. oscilla- 
tions in vacuum. The best-fit values of the oscillation parameters which were 
found from the analysis of the data 

(Am2)KL = 6.9 ■ lO-^eV^; (sin^ 2 ^)kl = 1- (9) 

are compatible with the solar LMA values. 

From common analysis of the solar and KamLAND data (assuming CPT) 
the following best-fit values of the parameters 

(Am2),„i+KL = 7.1 ■ lO-^eV^; (tan^ OU^kl = 0.41 (10) 

were obtained [3]. 

In from the global two-neutrino analysis of the solar and the Kam- 
LAND data the following 90% CL ranges of the parameters were found 

6.0-10-^ < (Am2),oi+KL < 8.7 -10" 0.25 < (sin^ 0)soi+kl < 0.37. (11) 

SK atmospheric neutrino evidence of neutrino oscillations was confirmed by 
the long baseline accelerator K2K experiment ^1]. The distance between 
neutrino source (KEK accelerator) and SK detector is about 250 km. Average 
energy is 1.3 Gev. In the K2K experiment 56 z/^ events were observed. 
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The expected number of the events is equal to 80.ll5;4. The best-fit values 
of the two oscillation parameters found from the two-neutrino analysis of the 
data 

(sin2 2^)K2K = 1; (Am')K2K = 2.8 10"^ eV^ 

are compatible with the atmospheric values Q. 

The negative results of the CHOOZ yiSj and Palo Verde JHl reactor 
experiments are very important for the neutrino mixing. In these experiments 
reactor Ve 's are detected via the observation of the classical process z/g + p — »• 

+ n. The distances between reactors and detectors in these experiments 
were about 1 km. No disappearance of z/g's were observed. In the CHOOZ 
experiment for the ratio of the total numbers of observed and expected events 
the following value 

= 1.01 ± 2.4% ± 2.7% 

exp 

was found. From the CHOOZ two-neutrino exclusion curve at Am^ = 2 • 
10~^eV^ (the SK best-fit value) the following upper bound 

(sin^2^)cz < 2 ■ 10"^ (12) 

can be obtained. 



3 Neutrino oscillations in the framework of 
three-neutrino mixing 

Neutrino oscillation data are analyzed under three basic assumptions (see, 
for example, [T7]): 

• Three fiavor neutrinos i^e, ( and antineutrinos z/g, ,z/^ ) exist 
in nature. 

• The Lagrangian of the interaction of the flavor neutrinos is given by 
the Standard Model. The lepton charged current and neutrino neutral 
current have the form 

= 2 UiLlah] jT = Yl ^lL7am. (13) 
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• Neutrino mixing 

3 

l^lL = Uli ViL (14) 

i=l 

takes place. Here z/j is the field on neutrino with mass irii and U is the 
unitary PMNS IT^ mixing matrix. 

The three-neutrino scheme fll4|) is the minimal scheme of neutrino mixing: 
the number of the massive neutrinos z/j is equal the the number of the flavor 
neutrinos vi (three). If the number of light neutrinos z/j is larger than three, 
sterile neutrinos must exist. An indication in favor of neutrino oscillations 
with the third relatively large neutrino mass-squared difference (~ 1 eV^) 
was obtained in the single accelerator LSND experiment [201 • '^^^ LSND 
result needs confirmation. It will be checked by the MiniBooNE experiment 
at the Fermilab [21] . 

Neutrino oscillation data are compatible with two types of neutrino mass 
spectra^. 

1. Normal mass spectrum (NS): 

mi < ms < mg; ^ml^ ~ Am,^„i_KL; ^"^32 - ^"^sk- 

2. Inverted mass spectrum (IS): 

m-i < nil < ^^21 — ^"^soI-kl! ^"^31 — ^^'^SK 

Here Amgj^ and Am'^^i_y^i^ are neutrino mass-squared differences, which 
are determined from the two-neutrino analysis of the atmospheric SK and 
the solar-KamLAND data (see Q. (IHl) . (ITIH) and (ITTll ). 

From neutrino oscillation data it follows that two independent neutrino 
mass-squared differences satisfy the inequality: 

Amli < Am^2- (15) 
In the framework of the three-neutrino mixing we have (see pTj ) 

4|f/,3p(l-|t/e3|') = (sin2 2^)cz. (16) 

^Ar7i^.j = mf, — mf. Notice that different notations for neutrino masses are used for NS 
and IS spectra 
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If we take into account solar data, from (|12|) and we obtain that 



U,.? < 5-10 



-2 



(17) 



The present-day picture of neutrino oscillations is determined by the in- 
equahties (fTH|) and (fTTjl (see, for example, |IZ|). In the leading approxi- 
mation neutrino oscillations in the experiments with L/i? ~ 10'^ (L is the 
source-detector distance in m (km), E is the neutrino energy in MeV (GeV)) 
are two-neutrino —>■ Vr oscillations. Solar neutrino transitions in matter 
(and antineutrino oscillations in reactor experiments with LjE 'r^ 10^) are 
^ii,T {pe ^ii.T) oscillations. Solar z/g and reactor v^. survival probabili- 
ties depends on two parameters fSTn^^ ^^^d sin^ Qyi (if CPT is assumed) and 
are given by the standard two-neutrino expressions in matter and in vacuum, 
respectively. 

The observation of neutrino oscillations means that the flavor lepton num- 
bers Lg, and are not conserved by a neutrino mass term. There exist 
two theoretical possibilities in this case 

1. The total lepton number L = + + is conserved. In this case 
fields of neutrinos with definite masses (x) are Dirac fields of neutrinos 
(L=l) and antineutrinos (L=-l). 

2. The total lepton number L is violated. In this case the fields of neu- 
trinos with definite masses i^i{x) are Majorana fields of neutrinos (with 
antineutrinos identical to neutrinos). The fields Uilx) satisfy the Ma- 
jorana condition 



where C is the charge conjugation matrix. 

After the discovery of the neutrino mixing the problem of the nature of 
massive neutrinos (Dirac or Majorana?) is one of the most fundamental. 
There is no doubt that the solution of this problem will have very important 
impact on our understanding of the origin of small neutrino masses and 
neutrino mixing. 

The investigation of the neutrino oscillations is an extremely sensitive 
method of the measurement of the very small neutrino mass-squared differ- 
ences [221 • This method does not allow, however, to reveal the nature of 
massive neutrinos |23j . 




(18) 



7 



In fact mixing matrices in the Majorana and Dirac cases are connected 
by the relation 

where Sik{l3) = e*^'' 6ik is a diagonal phase matrix. The matrix S{(3) does 
not give contribution to the transition probability 



Thus, we have 



P^^(z/„ ^ = P^(z/„ z/„0 (19) 



The nature of the massive neutrinos can be revealed only via the investigation 
of the processes in which the total lepton number L is not conserved. The 
most sensitive to small Majorana neutrino masses process is neutrinoless 
double p- decay of even-even nuclei (see reviews [2 ^ 125 } l26| l27 t 128 } l29 | IHH]): 

(A,Z) ^ (A,Z + 2) +e- +e-. 

In the case of the small Majorana neutrino masses the half-life of this 
process is given by the following general expression 

^ |meenMO'^(A,Z)|2G°'^(Eo,^). (20) 



Here |mee| is the effective Majorana mass 



imp 



(21) 



M'^'^{A, Z) is nuclear matrix element, which depends only on nuclear proper- 
ties, and G^'^{Eq, Z) is known phase-space factor {Eq is the energy release). 
In the next sections we will discuss 

1. Existing data and future experiments. 

2. Possible values of the effective Majorana mass |mee|. 

3. The problem of nuclear matrix elements. 
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4 Existing Oz//?/3-decay data and future exper- 
iments 



Neutrinoless double (3 decay is allowed for such even-even nuclei for which 
usual 13 decay is forbidden by the conservation of energy. There are several 
nuclei of this type: 

^'^Ge(2.039), ^^°Te(2.528), i=^'^Xe(2.480)/°° Mo(3.034)/^° Nd(3.367) 

and others. In the brackets the energy release Eq in MeV is given. This 
is an important characteristic of the Oz//3/3-decay: the decay probability is 
proportional to E^. 

The results of many experiments on the search for Qv(3(3 -decay are avail- 
able at present (see jSHESESI)- The most stringent lower bounds on the half- 
life of p(3 /3-decay was reached in the Heidelberg-Moscow and CUORI- 
CINO jHSl experiments. 

The detector (and source) of the Heidelberg-Moscow experiment ^341 con- 
sists of 5 crystals of 86% enreached ''^Ge of the total mass about 11kg . For 
the half-life the lower bound 

Tf/s > 1-9 ■ 10^^ y (90% CL) (22) 

has been found. Taking into account different calculations of the nuclear 
matrix element, from (|22j) for the effective Majorana mass |mee| the following 
upper bounds 

\mee\ < (0.3-1.2) eV. (23) 

were obtained. 

In the cryogenic experiment CUORICINO Te02 crystals with a total 
mass 40.7 kg are employed. For the half-life of ^^''Te in this experiment the 
following lower bound 

Ti% > 7.5 ■ 10^^ years (24) 

was reached recently. From ()24|) for the effective Majorana mass the upper 
bounds 

l^eel < (0.3-1.7) eV (25) 

were obtained. 

Many projects of new experiments on the search for the neutrinoless dou- 
ble /9-decay of different nuclei are under research and development at present 
(see j21ElEn| )• 
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The main goal of the future experiments is to reach the sensitivity |mee| — 
a fewlO~^eV. This goal can be accomplished by detectors with mass about 
1 ton or more, which have a good energy resolution, low background and an 
efficient signature for Ou/SP events. 

The experiment CUORE [33] will be a continuation of the CUORICINO 
experiment. Cryogenic detector will consist of 1000 Te02 crystals operated 
at a temperature 10 mK. The total mass of the detector will be about 800 
kg. The expected resolution at Eq = 2.528 MeV is 5 keV. For the half-life of 
^^'^Te the value 

r°/2 - 9.5- 10^^ years (26) 
is envisaged. This corresponds to the sensitivity 

\mee\ ^ (2-5.2) -10-2 eV (27) 

In the EXO experiment up to 10 tons of 60-80 % enreached ^^^Xe are 
planned to use. An important feature of this experiment is a laser tagging 
of Ba^ ions, produced in the recombination of Ba^"*" ions from the decay 
^^^Xe — i>^'^6 Ba^"*" -|- e~ -|- e~. The detection of Ba"^ ions will provide large 
background reduction. The value 

r°/2 - 1 ■ 10^® years (28) 

is expected. It corresponds to the sensitivity 

|mee| ~ (1.3-3.7) ■ 10"^ eV (29) 

The GENIUS experiment |HZ1 will be a development of the Heidelberg- 
Moscow experiment. About 1 ton of 86 % enreached ^^Ge will be embadded 
in a large liquid nitrogen cryostat . The liquid nitrogen will provide effective 
shielding from the external background. For the half-life a value 

T^/a - 1 ■ lO^Vears (30) 

is expected. This value corresponds to the sensitivity 

|mee| ~ (1.3-5.0) ■ 10"^ eV (31) 

In the MAJORANA experiment which will be the continuation of 
the IGEX experiment jHH], about 500 kg of 86 % enreached ^^Ge will be used. 
The main background is expected from the decay ^^Ge -^^^ Ga -|- e"*" -|- u^. 
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It will be suppressed by the segmentation of the detector and effective pulse 
shape analysis of the' signal. In the MAJORANA experiment the value 



T 



1/2 



4 ■ 10 years 



(32) 



is expected. It corresponds to the sensitivity 



m, 



'ee 



(2.1-7.0) ■ 10-2 eV. 



(33) 



5 Effective Majorana mass 

The effective Majorana mass is determined by the absolute values of the 
neutrino masses, mixing angles and CP phases. Let us discuss first these 
three ingredients. 

1. Neutrino masses 

From neutrino oscillation data only neutrino mass-squared differences 
can be inferred. In the case of the NS neutrino mass spectrum the 
neutrino masses are given by the relations 



From existing data only upper bounds of the lightest mass mmin(mi or ms) 
can be obtained. From the data of the Troitsk [lU] and Mainz |39| tri- 
tium /3-decay experiments the following bounds were found 



More stringent bound 

mmin < 0.6 eV 

was obtained from the analysis jH] of the data of WMAP and Sloan 
Digital Sky Survey Collaborations. 




For the IS neutrino mass spectrum we have 




^min < 2.05 (2.3) eV 



(34) 
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2. Mixing angles. 

The elements Uei (i=l,2,3) are given by 

^ei I ^ei \ ^ 7 

where ctj are Majorana CP-phases. If we introduce the angle ^13 in 
such a way that 

|[/e3p = sin2^i3 
from the unitarity relation X]^=i l^^iP = 1 we have 

|[/eiP = cos^6'i3 cos^6'i2; |?7e2p = cos^ 6*13 siu^ 6*12. 

The parameter sin^ 9i2 can be determined from the results of the solar 
and KamLAND experiments. From the analysis of the existing data 
the best fit value (fTIH) and a range (fTTj) were found. 

Only upper bound of the parameter sin^ ^13 is known today. From the 
exclusion curve obtained in the CHOOZ experiment the bound ()17p 
was obtained. 

3. Majorana phases 

Majorana phases are unknown. In the case of the CP invariance in 
the lepton sector the elements of the mixing matrix satisfies the relation 

m 

Ue^ = U:,V^, (35) 

where rji = i pi = ±1 is the CP-parity of the Majorana neutrino Ui. 
Thus, in the case of the CP invariance 

and the effective Majorana mass takes the form 

3 

Imeel = l^lUeil"^ Pi TUil (36) 
i=l 

From this equation it follows that in the case of the different CP par- 
ities of Majorana neutrinos with definite masses cancellations of their 
contributions to the effective Majorana mass can take place. In the 
general case of the CP violation \mee\ depends on two Majorana phase 
differences. 
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The value of the effective Majorana mass |mee| strongly depends on the 
pattern of the neutrino mass spectrum and lightest neutrino mass PH" . Three 
types of the neutrino mass spectrum are usually considered. 

• Neutrino mass hierarchy ^ 



mi <S 1712 <C m^. (37) 

In the case of the neutrino mass hierarchy 1712 and rris are determined 
by the "solar-KamLAND" and "atmospheric" mass-squared differences 



7712 — y Amli, 7n^ ~ y Aml2 (38) 

and the lightest neutrino mass mi is small: mi ^ ^yAm^ ~ 8.4 • 
lO-^eV. 

Neglecting the contribution of mi to the effective Majorana neutrino 
mass, we have 



cos^ ^13 sin^ 612 \/ Am|i + e'^"''^ sin^ ^13 J Am^ 



''32 



(39) 



where 032 = 03 — 0:2. It follows from ()39|) that \mee\ is small: the first 
term is small because of the smallness of ^/Am^; the contribution of 
the "large mass" ^/Am^ is suppressed by the smallness of the param- 
eter sin^ 6'i3. From Q (fTT|) (fTTj) and we obtain the following upper 
bound (90 % CL): 



|mee| < 6.0- 10~W (40) 

The bound ()40|1 is significantly smaller than the sensitivity of the future 
experiments on the search for Oz//3/3-decay. The observation of the 
neutrinoless double /3-decay in the experiments of the next generation 
would mean that the neutrino masses do not follow the hierarchy ()37|). 

• Inverted hierarchy of neutrino masses 

m^ <^ mi < m2- (41) 

■'Notice that masses of quarks (up and down) and charged leptons follow the hierarchy 
of the type (p?7jl . 
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In the case of the inverted mass hierarchy mi and m2 are determined 
by "atmospheric" mass-squared difference 



mi ~ m2 — y |Am3i| (42) 

and the lightest neutrino mass ma is small: ^ y^|Am|J ^ 4.5 • 
10-2eV. 

Neglecting small contributions of and sin^ 613, for the effective Ma- 
jorana mass we have 



^^3il I Yl - V (l-sin'2^12 sin2a2i)^ (43) 



1=1.2 



where 0^21 = 02—0:1. In the case of the CP invariance 0:21 = f (p2— Pi) = 
±- 

u, _i_ 2 . 

From Eq. ()43p for the effective Majorana mass we have the following 
range 



cos2 0i2 Y |Am|i| < |mee| < y\Aml^\, (44) 

where the upper and lower bounds correspond to the cases of the CP 
conservation: the upper bound corresponds to the case of equal CP 
parities of z/i and 1^2 and the lower bound corresponds to the case of 
opposite CP parities. 

From ()44j) and (jllj) at 90% CL for the effective Majorana mass 
we have the range 



m 



ee 



0.26^|Am|i| < |mee| < yl^^iil (45) 

Thus, in the case of the inverted hierarchy of neutrino masses the scale 
of the effective Majorana mass is determined by A/|Am|7[. From 
and (jnj) we have 

0.9-10"W< |mee| <5.5-10"W (46) 

The values of the effective Majorana mass in the range (j^6|) can be 
reached in Oi//5/?-decay experiments of the next generation. Let us stress 
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that in the case of the inverted hierarchy of the neutrino masses the 
value of |mee| can not be smaller than ~ 1 ■ 10~^eV. This is connected 
with the fact that the value sin^ 2 9i2 = 1 (maximal 1-2 mixing) is 
excluded at 5.4 a level by the solar and KanLAND data [3]. 

From Eq. (j43j) we find 



sin^a^i^ (l-^^l . (47) 

If the problem of the nuclear matrix elements will be solved (see the 
next section for discussion), the experiments on the measurement of 
the half-life of Oi>f3(3-decaj could allow to obtain an information on the 
the value of the CP parameter sin^ 

Practically degenerate neutrino masses 
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If the lightest neutrino mass is much larger than ^^/Arn^ ^ 4.5-10 ^eV 
in this case neutrino masses are practically degenerate: 

mi ~ m2 — m^. (48) 

For the degenerate neutrino masses the effective Majorana mass 

3 

\mee\ — riiQ |^f/g^-|~mo (1 — sin^ 2 6'i2 sin^ a) 2 , (49) 

i=l 

depends on two parameters: sin^ a {a is the Majorana CP phase dif- 
ference) and a common mass rriQ. From ()49|1 and pijl we have the 
range 

0.26 mo < |mee| < mo (50) 

Thus, in the case of the practically degenerate Majorana neutrino mass 
spectrum the scale of the effective Majorana mass |mee| is determined 
by (unknown) common mass mo. 

The mass mo can be determined from the data of the experiments on the 
measurement of the high-energy part of the /9-decay spectrum of tritium 
and from cosmological data. In the tritium experiment KATRIN p9| 
1^ . now at preparation, the sensitivity mo — 0.2 eV is expected. 

If in the future Oi//?/?-decay experiments it will be found that the value 
of the effective Majorana mass is significantly larger than a/ Amgg — 
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4.5 ■ 10 ^eV from Oi//3/?-decay data an information about the value of 
the common mass tuq can be inferred: 

\mee\ <mo < ^"1^^^ (51) 
cos 2 0i2 

Using the existing data (see (fTTj) ) we have 

\mee\ < Tno < 3.8 Imeel 



In conclusion we would like to emphasize that the measurement of 
the effective Majorana mass |mee| could allow to obtain an important 
information on the pattern of the Majorana neutrino mass spectrum. 



6 The problem of the nuclear matrix elements 

Neutrinoless double /?-decay is a second order in the Fermi constant Gp 
process with a virtual neutrino. For small neutrino masses (much smaller 
than the bounding energy of nucleons in nuclei) the matrix element of the 
OuPP -decay is factorized in the form of a product of the effective Majorana 
mass I rriee \ , which depends on neutrino masses and elements U^^ of the mixing 
matrix, and nuclear matrix element (NME), which is determined only by 
the strong interaction. The NME is a matrix element of the chronological 
product of the two CC hadronic currents and the neutrino propagator. It can 
not be connected with other observables. In the calculation of NME many 
intermediate states must be taken into account. 

Two basic methods of the calculations of NME are used : quasiparti- 
cle random phase approximation (QRPA) and nuclear shell model (NSM). 
Many calculations of NME of different nuclei, based on these approximate 
approaches, exist in literature (see [211 123 1221) • The results of different cal- 
culations of NME differ by about factor three or more. For example, if we 
assume that \mee\ = 5 ■ 10"^ eV from different calculations of NME for the 
half-life of the Oz//3/3-decay of ^^Ge the values in the range 

6.8 ■ lO^V < T°;;2(^^Ge) < 70.8 ■ lO^V 
can be obtained [221 • 
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Recently a progress in the calculation of NME in the framework of QRPA 
have been achieved HB]. The nuclear matrix elements of Qvj3j3- decay of ''^Ge, 
^°°Mo, ^^°Te and ^^^Xe were calculated with the values of the parameter of 
the effective particle-particle interaction Qpp determined from the measured 
life-time of the 2i> j3l3- decay. It was shown that the values of nuclear matrix 
elements of the Oz//5/5-decay of these nuclei are stable under the change of 
the nuclear potential and the number of single particle states used as a ba- 
sis. Moreover the matrix elements calculated in different QRPA models are 
practically the same (differ not more than 10%.) 



7 Possible test of the models of NME calcu- 
lations 

Taking into account all uncertainties connected with the calculations of the 
nuclear matrix elements of the Oz//3/3- decay, it will be very important to 
find a possibility to test NME calculations. We will discuss here a possible 
test, based on the factorization property of the matrix elements of the Oz//3/3- 
decay (see^j). 

The proposed test can be realized if neutrinoless double /?-decay of several 
nuclei Ai, Zi (i=l,2,...) is observed. Using a model M of the calculation 
of the nuclear matrix elements, from Eq. (j2Uj) the value of the parameter 
\mee\% z,(^) can be determined. The model M is compatible with the data 
if the relations 

\mee\X,zAM) - |mee|^„^,(M) = ... (52) 

are satisfied. 

From (^01) it follows that for any nuclei (A,Z) the product 



\m,,\lziM)\M'^iA,Z)\ 



M 



does not depend on the model M. Thus, for two different models Mi and M2 
we have 

\mee\lziM2) = |mee|iz(Mi)77*^^'*^HA,Z), (53) 



where 



^ ^^'^)-|Mo^(A,z)ii,; 
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In the Table I we present the values of the coefficient ri{A, Z) for the case of 
the matrix elements calculated in (NSM) and in <A6, (the latest QPRA 
calculations). 

If ri{A, Z) depends on (A, Z) and one model is compatible with the data 
the other model in principle can be excluded. However, as it is seen from the 
Table I from the observation of the Oz//?/?- decay of ^^^Xe and ^^°Te it will be 
difficult to distinguish models |1H] and jlH]: the difference between 7][^^^X.e) 
and ?7(^^°Te) is about 20% . It will be more easier to distinguish models jlHl 
and if Oz//3/3- decay of ^^Ge and ^^°Te is observed. 

Taking into account the existence of many models of the calculations of 
the nuclear matrix elements of the Oz//3/3- decay we can conclude that the 
observation of neutrinoless double /3-decay of three (or more) nuclei would 
be an important tool in the solution of the problem of NME. 

Table I 

The parameter ■q^SM;QRPAi^j^^ determined by Eq. (j33)), for nuclear 
matrix elements of the Oz//3/3- decay, calculated in Ref.^j (NSM) and in 

Ref.jini (QRPA) 



Nucleus 


^NSM;QRPA 


^^Ge 


3.1 




2.1 


i36Xe 


2.5 



8 Conclusion 

After the discovery of neutrino masses and neutrino mixing the problem of the 
nature of neutrinos with definite masses z/j is one of the most fundamental. 
The establishment of the nature of z/j will have a profound impact on the 
understanding of the mechanism of the generation of the neutrino masses 
and mixing. 

The most sensitive to the small Majorana neutrino masses process is neu- 
trinoless double /?-decay of even-even nuclei. From today's data to following 
bounds for the effective Majorana mass can be inferred 

\mee\ < (0.3 - 1.2) eV 

New experiments on the search for Oz//3/?-decay of ^^°Te, ^^Ge, ^^^Xe, ^"^Mo 
and other nuclei are in preparation at present. In these experiments the 
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sensitivity 

\mee\ — a fewlO"^ eV 

is envisaged. 

The data of neutrino oscillation experiments allow to predict ranges of 
possible values of the effective Majorana mass for different patterns of the 
neutrino mass spectra. In order to obtain information on the neutrino mass 
spectrum it is important not only to observe Oz//3/3-decay but also to deter- 
mine the value of the effective Majorana mass \mee\- 

From the measured half-life of Oz//3/?-decay only the product of the effec- 
tive Majorana mass and nuclear matrix element can be obtained. Existing 
calculations of the nuclear matrix elements of the Oz//5/5-decay differ by about 
a factor of three or more. The improvement of the calculations of the nuclear 
matrix elements is a real theoretical challenge. We have discussed here a pos- 
sible method which could allow to test models of calculation of the nuclear 
matrix elements of the Oz//?/?-decay. The method is based on the factorization 
property of the matrix element of 0z//3/3-decay and require observation of the 
Oz//3/3-decay of several nuclei. 

I would like to thank S.T. Petcov for usefull discussions. I acknowledge 
the support of the Italien Program "Rientro dei cervelli" . 
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